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DETERMINING THE FIRST ORDER PERTURBATION OF A 
POLYHARMONIC OPERATOR ON ADMISSIBLE MANIFOLDS 


YERNAT M. ASSYLBEKOV AND YANG YANG 


Abstract. We consider the inverse boundary value problem for the first or¬ 
der perturbation of the polyharmonic operator £ g .x.q- with X being a V! 71 ’ 00 
vector field and q being an L°° function on compact Riemannian manifolds 
with boundary which are conformally embedded in a product of the Euclidean 
line and a simple manifold. We show that the knowledge of the Dirichlet-to- 
Neumann determines X and q uniquely. The method is based on the construc¬ 
tion of complex geometrical optics solutions using the Carleman estimate for 
the Laplace-Beltrami operator due to Dos Santos Ferreira, Kenig, Salo and 
Uhlmann. Notice that the corresponding uniqueness result does not hold for 
the first order perturbation of the Laplace-Beltrami operator. 


1. Introduction 


Let (M, g) be a compact oriented Riemannian smooth manifold with boundary. 
Throughout this paper, the word “smooth” will be used as the synonym of “(7°°”. 
Let A g be the Laplace-Beltrami operator associated to the metric g which is given 
in local coordinates by 



where as usual is the matrix inverse of (gjk), and \g\ = det (gjk)- If F denotes 
a function or distribution space (C k , L p , H k , V. etc.), then we will denote by 
F(M,TM ) the corresponding space of vector fields on M. 

Let X £ W 1,ao (M,TM) and q £ L°°(M). Consider the polyharmonic operator 
(—A g ) m , to > 1, with the first order perturbation induced by X and q 


£g,x,q — (-A g) m + X + q 
The operator C g ,x,q equipped with the domain 

F>(£ g ,x, q ) = {«£ H 2m (M) : yit = 0} = H 2m (M) n H™{M) 


is an unbounded closed operator on L 2 (M ) with purely discrete spectrum; see [8]. 
Here and in what follows, 


7 U := (u\dM, A g u|aM, • • •, A^ 1 1 u\g M ) 


is the Dirichlet trace of u, and H S (M ) is the standard Sobolev space on M, s £ R. 
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We make the assumption that 0 is not a Dirichlet eigenvalue of £ g ,x, q in M. Under 
this assumption, for any / = (/o, - - -, fm-i) £ U m {dM) := rij’L’o 1 H 2rn ~ 2 ^ x / 2 {dM ), 
the Dirichlet problem 

£g,x, q u = 0 in M, 
ju = f in d M, 

has a unique solution u £ Let v be an outer unit normal to dM. Intro¬ 

ducing the Neumann trace operator 7 by 

m— 1 

7 : H 2m {M) H 2m ~ 2: >- 3/2 {dM), 

j =0 

7 U = ( d v u\dM , d„A g u\g M , ■ ■ ■, d v A™ l u\g M ), 

we define the Dirichlet-to-Neumann map N Qi x,q by 

771—1 

N g , x ,q : n m (dM) -»• J] H 2m - 2 i~ 3 ' 2 (dM), N g , x , q (f) = 7 u, 

3 =0 

where u £ H 2m (M) is the unique solution to the boundary value problem (1). Let 
us also introduce the set of the Cauchy data for the operator C g .x. q 

C g ,x,q = {{'yu, 7 yu) : u £ H 2m (M ), £ g ,x, q u = 0}. 

When 0 is not a Dirichlet eigenvalue of C g< x,q in M, the set C gt x, q is the graph of 
the Dirichlet-to-Neumann map N gyX , q - 


The inverse problem we are concerned in this paper is to recover the vector field X 
and the function q from the knowledge of the Dirichlet-to-Neumann map N g ,x, q on 
the boundary dM. 

When rri = 1, the Dirichlet-to-Neumann map N 9jX , q is invariant under gauge trans¬ 
formations in the following sense. Let ip be a C 2 [M ) such that ip\dM = 0 and 
d„ip\dM = 0. Then 


4 J0. g ,x, q ^ = C 


g,x,q’ 


^g,X,q ~ Ng,X,q, 


where 

X = X + 2V'0, q = q+{X,Vif>) g + \Vil>\ 2 g -iA g i(>. 

Therefore, we may hope to recover X and q from boundary measurements only 
modulo the above gauge transformations. 

In the Euclidean setting, this inverse boundary value problem has been extensively 
studied, usually in the context of magnetic Schrodinger operators [15, 18, 19, 22, 24]. 
In the case of Riemannian manifolds, this was proved in [3] for the special class of 
so-called admissible manifolds. 


Let us now introduce admissible manifolds. For this we need the notion of sim¬ 
ple manifolds [21]. The notion of simplicity arises naturally in the context of the 
boundary rigidity problem [17]. 
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Definition 1.1. A compact Riemannian manifold ( M,g ) with boundary is said to 
be simple if the boundary dM is strictly convex, and for any point x G M the 
exponential map exp^ is a diffeomorphism from its maximal domain in T X M onto 


M. 


Definition 1.2. A compact Riemannian manifold ( M,g ) with boundary of dimen¬ 
sion n > 3, is said to be admissible if it is conformal to a submanifold with boundary 
of 1 x (Mo, go) where (Mo, go) is a simple (n — l)-dimensional manifold. 

Examples of admissible manifolds include the following: 

1. Bounded domains in Euclidean space, in the sphere minus a point, or in hy¬ 
perbolic space. In the last two cases, the manifold is conformal to a domain in 
Euclidean space via stereographic projection. 

2. More generally, any domain in a locally conformally flat manifold is admissible, 
provided that the domain is appropriately small. Such manifolds include locally 
symmetric 3-dimensional spaces, which have parallel curvature tensor so their Cot¬ 
ton tensor vanishes (see the [3, Appendix B]). 

3. Any bounded domain M in R™, endowed with a metric which in some coordinates 
has the form 



with c > 0 and go simple, is admissible. 

4. The class of admissible metrics is stable under C 2 -small perturbations of go- 


It was shown in [12] that, in the Euclidean case, the obstruction to uniqueness 
coming from the gauge equivalence when m = 1 can be eliminated by considering 
operators of higher order. The purpose of this paper is to extend this result for the 
case of admissible manifolds. 

Our main result is as follows. 

Theorem 1.3. Let (M,g) be admissible, and letm> 2 be an integer. Suppose that 
X 1 ,X 2 G lf 1 , 00 (KxMo,r(lxMo))n^'(M, TM) and qi ,q 2 G L°°(M) are such that 
0 is not a Dirichlet eigenvalue of Cg^x-i^ and C g ,x 2 ,q 2 * n M. If Ng^Xi^ = Ng,x 2 ,q 2 > 
then X\ = X 2 and q\ = q 2 . 

The key ingredient in the proof of Theorem 1.3 is the construction of complex geo¬ 
metric optics solutions for the operator C g ,x,q with X being a IT 1 ’ 00 vector field 
and q an L°°(M) function. For this, we use the method of Carleman estimates 
which is based on the corresponding Carleman estimate for the Laplacian due to 
Dos Santos Ferreira, Kenig, Salo and Uhlmann [3]. 

In Theorem 1.3, the condition that X\ = X 2 = 0 on dM is needed to extend the 
vector fields X\ and X 2 to a slightly larger simple manifold than M while preserving 
the W 1 ' 00 regularities. When more regularities on Xj and qj (j = 1, 2) are available, 
we can show a boundary determination result for the vector fields and thus drop 
such an assumption. This is the following theorem. 
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Theorem 1.4. Let (M,g) be admissible, and let m > 2 be an integer. Suppose 
that Xi,X 2 £ C°°(M,TM ) and qi,q 2 £ C°°(M ) are such that 0 is not a Dirichlet 
eigenvalue of C g ^x!, qi and C g ,x 2 ,q 2 in Af. If N 3i Xi,qi — ^g,x 2 ,q 2! then X\ = X 2 
and qi = < 72 - 

Let 7r : H. x M$ —» M$ be the canonical projection 7r(xi,a/) = x'. It is interesting 
to notice that the boundary determination becomes unnecessary if go) is a 

simple (n — l)-dimensional manifold and dM is connected. 

Theorem 1.5. Let ( M,g ) be admissible, and letm> 2 be an integer. Suppose that 
Xi , X 2 £ W 1 ’°°(M,TM) and < 71 , <72 £ L°°(M) are such that 0 is not a Dirichlet 
eigenvalue of L g Xx,q\ an d £g,x 2 ,q 2 In M. Suppose further that (ir(M),g 0 ) is a 
simple (n — 1) -dimensional manifold and dM is connected. If N g ^Xi, qi = Ng,x 2 ,q 2 > 
then X\ = X 2 and qi = 92 - 

In the case of Euclidean space, the recovery of a zeroth order perturbation of the 
biharmonic operator, that is when m = 2, has been studied by Isakov [11], where 
a uniqueness result was obtained, similarly to the case of the Schrodinger operator. 
The recovery of a first order perturbation of the biharmonic operator from partial 
data was studied in [13] in a bounded domain, and in [25] in an infinite slab. Higher 
order operators occur in the areas of physics and geometry such as the study of the 
Kirchhoff plate equation in the theory of elasticity, and the study of the Paneitz- 
Branson operator in conformal geometry; for more details see [7]. 

Finally, we would like to remark that the problem considered in this paper can be 
viewed as generalization of the Calderon’s inverse conductivity problem [1], known 
also as electrical impedance tomography. In the fundamental paper by Sylvester and 
Uhlmann [23] it was shown that C 2 conductivities can be uniquely determined from 
boundary measurements. A corresponding result was proved by Dos Santos Ferreira, 
Kenig, Salo and Uhlmann [3] in the setting of admissible geometries. 

The structure of the paper is as follows. In Section 2 a Carleman estimate is 
derived for polyharmonic operators based on a similar estimate for the Laplace- 
Beltrami operator. Section 3 is devoted to the construction of complex geometric 
optics solutions for the perturbed polyharmonic operator C g ^x,q with X being a 
W 1 ’ 00 vector field and q £ L°°(M). Then the proof of Theorem 1.3 is given in 
Section 4. Attenuated ray transform is the subject of Section 5. In Section 6, we 
show that the Dirichlet-to-Neumann map determines X on the boundary, this leads 
to the proof of Theorem 1.4. Finally, the proof of Theorem 1.5 is given in Section 
7. 


2. Carleman estimates for polyharmonic operators 

Let (M, g) be a Riemannian manifold with boundary. In this section, following [3, 
14] , we shall use the method of Carleman estimates to construct complex geometric 
optics solutions for the equation C g ^x, q n = 0 in M, with X being a W 1,ao vector 
field on M and q £ L°°(M). 

We start by recalling the definition of the Carleman weight for the semiclassical 
Laplace-Beltrami operator —h 2 A g . Let U be an open manifold without boundart 
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such that M CC U and let ip C C°°(U, R). Consider the conjugated operator 

P v = e ^ h {-h 2 Ag)e~^ h . 

Following [3, 14], we say that ip is a limiting Carleman weight for —h 2 A g in U , if it 
has non-vanishing differential, and if it satisfies the Poisson bracket condition 

{p^,Pip}{x,0 = 0 when p v (x,£) = 0, (x,£)eT*M, 

where p v is the semiclassical principal symbol of P v . 

First we shall derive a Carleman estimate for the semiclassical polyharmonic oper¬ 
ator (—h 2 A g ) m , where h > 0 is a small parameter, by iterating the corresponding 
Carleman estimate for the semiclassical Laplace-Beltrami operator —h 2 A g , which 
we now proceed to recall the following [3, 14]. 

We use the notation dVol g for the volume form of (M,g). For any two functions 
u,v on M, define an inner product 

(u\v) := / u(x)v(x) dVo\ g (x), 

Jm 

and the corresponding norm will be denoted by || • ||l 2 (m)- We also write for short 

l|Vu|| L2(M) = |||Vu||| L2(M) = (jjVu(x)\ 2 g dVo\ g (x)^j ' . 

We assume that (M, g) is embedded in a compact manifold (N, g) without boundary, 
and tp is a limiting Carleman weight on ( [U,g ), where U is an open submanifold of 
N such that M CC U. By semiclassical spectral theorem one can define for s £ M 
the semiclassical Bessel potentials J s = (1 — h 2 A g ) s / 2 . One has J s J l = J s+t , 
and Bessel potentials commute with any function of — A g . Define for s £ R the 
semiclassical Sobolev space associated to the norm 

= \\J s u\\mN)- 

Our starting point is the following Carleman estimate for the semiclassical Laplace- 
Beltrami operator —h 2 A g which is due to Dos Santos Ferreira, Kenig, Salo and 
Uhlmann [3, Lemma 4.3]. In what follows, A < B means that A < CB where 
C > 0 is a constant independent of h and A, B. 

Proposition 2.1. Let ( U,g) be an open Riemannian manifold and ( M,g ) be a 
smooth compact Riemannian submanifold with boundary such that M CC U. Let 
p be a limiting Carleman weight on ( U,g). Then for all h > 0 small enough and 
s£l, we have 

^hll^+biv) £ \\ eV/h (-h 2 A g ) e -^ hu \\ HLi{N) 

for all u C Cq°(M). 

Next we shall derive a Carleman estimate for the operator C g .x,q with X being a 
W 1,00 vector field on M and q G L°°(M). To that end we shall use Proposition 2.1 
with s = — 1. We have the following result. 
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Proposition 2.2. Let ( U,g ) be an open Riemannian manifold and (M,g ) be a 
smooth compact Riemannian submanifold with boundary such that M CC U. Let ip 
be a limiting Carleman weight on ( U,g ). Suppose that X is a W 1,0 ° vector field on 
M and q £ L°°(M). Then for all h > 0 small enough, we have 


IMU’W < j^\\ eV/h {h 2m C g ,x,q)e-*/ h u\\ H j {NV 

(2) 

for all u £ 


Proof. Iterating the Carleman estimate in Proposition 2.2 m times, 
the following Carleman estimate for the polyharmonic operator, 

m > 2, we get 

h m M H s+m (N) < \\e*/ h (-h 2 A g ) m e-*/ h u\\ HUN) , 


for all u £ Cq°(L l), s £ R. and h > 0 small enough. We shall use this estimate with 
s = —1: 

nMliC-w iS \\ e * ,h (-h 2 ^ 9 ) me -* /h u\\H- l{N y (3) 

for all u £ C^°(N) and h > 0 small enough. Since we are dealing with first order 
perturbations of the polyharmonic operator and m > 2, the following weakened 
version of (3) will be sufficient for our purposes 

nMU’W ;$ \\e*' h (-h a A g re-*/ h u\\ H £ W , 

(4) 

for all u £ and h > 0 small enough. 

It is easy to see that 


\\ e ‘ p / h h 2rn qe~‘ p / h u\\ L 2 ^ N) < h 2m \\q\\ L - {M) \\u\\ HLi{N) . 

(5) 


Note that e‘ p ^ h h 2m X(e v ^ h u) = —h 2m 1 (X, Xip) g u + h 2m (X, Xu) g . Therefore, 
since m > 2 

Vip) g u\\ L 2 (N) < h 2m ~ 1 \\(X, V<p) g \\L°°(M)\\u\\L 2 (N) 

< h m \\(X,Xp,) g \\ Lao{M) \\u\\ HLi{N) 

and 

\\h 2m (X,hVu) g || L 2 (JV) < /i 2m_1 ||X|| i a 0 (M) IMIet^jv) < h m \\X\\ L ~ {M) \\u\\ HLiiN) . 
imply 

\\e*/ h h 2m X(e-^ h u)\\ LHN) < h m \\u\\ HUN) . 

Combining this together with estimates (4) and (5), we get the result. □ 

Set 

:=e^ h {h 2m C gtXtq )e-^ h . 

Then we have 

(£ v u,u)n = (u, £*u)n, u,veC£°(n), 

where £* = e~‘ p ^ h (h 2 C g ^x.-div g x+q)e ip ^ h is the formal adjoint of C<p, and (■, -)m 
is the distribution duality on M. The estimate in Proposition 2.2 holds for £*, 
since —ip is a limiting Carleman weight as well. 

To construct the complex geometric optics solutions for the operator C g ^x,q, we need 
to convert the Carleman estimate (2) for £* into the following solvability result. 
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The proof is essentially well-known, and we include it here for the convenience of 
the reader. We shall use the following notation for the semiclassical Sobolev norm 
on M 

IMIh s i c i(m) = IMIl 2 (m) + IIW U ||| 2(M ). 

Proposition 2.3. Let X be a W 1 ’ 00 vector field on M and q £ L°°(M) and assume 
that m > 2. If h > 0 is small enough, then for any v £ L 2 (M) there is a solution 
u £ H 1 (M) of the equation 

e V>Ih h 2m Cg ^ qe ->p/h u = y 

satisfying 

C 

Proof. Let v £ and let us consider the following complex linear functional, 

C, C;w^(w,v) M . 

By the Carleman estimate (2) for £*, the map L is well-defined. Let w £ Cq°(M). 
Then we have 

\L{C* v w)\ = \(w,v) M \ < |MU 2 (AT)IMIl 2 (M) 

By the Hahn-Banach theorem, we may extend L to a linear continuous functional L 
on L 2 (N), without increasing its norm. By the Riesz representation theorem, there 
exists u £ H 1 { K") such that for all £ H~ 1 ( R n ), 

= (ip,u) K n, and |M|tfi ci(R n) < —|H|L 2 ( n } . 

Let us now show that C^u = v in fl. To that end, let w £ Cq°(LI). Then 

(C v u,w) n = ( u , = L(C* V w) = (w,v)q = (v,w)n- 

The proof is complete. □ 

3. Complex geometric optics solutions 

Let ip be a limiting Carleman weight in an admissible manifold (M, g). We will 
construct solutions to C 9t x,qU = 0 in M of the form 

u = e -(p+#)/fe( a + r ), (6) 

where a is an amplitude, r is a correction term which is small when h > 0 is small, 
and if is a real valued phase. 

Set p = (p + iif for the complex valued phase. Consider the conjugated operator 
P p = e p / h h 2rn C gt x, q e- p/h , which has the following expression 

Pp = {-h 2 Ag - |Vp|g + hX g p + 2 hXp) m + h 2m X - h 2 ™- 1 (-■Vp, X) g + h 2m q. 
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Here and in what follows, the norm | • \ 2 and the inner product (•, •) g are extended 
to complex valued tangent vectors by 

(C,v)g = (Re(,Rer]) g -{lm(,lmri) g +i{(ReC,Iinri) g +{IrnC,Reri) g ), \(\ 2 g = (C,C )g- 

Since to > 2, in order to get 

e vl h h*mC gtX , q ( e -*/h a ) = 0(h m+1 ), 

in L 2 (M), we should choose p satisfying the following eikonal equation 

\V P \ 2 g = 0 in M, (7) 

and choose a C C°°(M ) satisfying the following transport equation 

(2Vp + A g p) m a = 0 in M. (8) 

Recall that ( M,g ) is conformally embedded in R x (Mo,go), where (Mo,go) is some 
simple (n — l)-dimensional manifold. If necessary, we replace Mo with a slightly 
larger simple manifold. Therefore, we can and shall assume that for some simple 
(D,g 0 ) CC (M) nt ,g 0 ) one has 

(M,g) CC (R x D int ,g) C (R x M‘ nt ,g). (9) 


Note that R x Mq has global coordinate chart in which the metric g has the following 
form 

9(x)=cm() (10) 

where c > 0 and go is simple. A natural choice of the limiting Carleman weight is 
<p(x) = X\. Then the equation (7) for the complex valued phase p becomes 

IW>| 2 = -, d xi il> = 0. 

c 

This equation will be solved using special coordinates on (M, g). This is based on the 
so-called polar coordinates on the transversal simple manifold (Mo,go)- Let ui G D 
be such that (x\,uj) £ M for all x\. Points of M have the form x = (x\ , r, 6) where 
(r, 9) are polar normal coordinates in ( D,go ) with center oj. That is, x' = exp ®(r9) 
where r > 0 and 9 £ S n ~ 2 . In terms of these coordinates the metric g has the form 

A ° 0 \ 

g(xi,r,8) = c(x\,r,9) I 0 1 0 , 

\0 0 m(r, 9)) 

where to is a smooth positive definite matrix. 

We solve (7) by simply taking ip( x ) = = r. Thus, the complex valued phase 

has the form p = X\ + ir and its gradient is Vp = -d, where 


Next, we solve transport equation (8). In the coordinates (x\,r, 9) equation (8) 
becomes 


bogM 


a = 0. 
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Consider a as the function having the following form 

a = \g\~ 1/4 c 1/2 a o (xi,r,0)b(6) 

where b is smooth and ao is such that dao = ca\ for some a\ satisfying da\ = 0. 
Note that (6) will be a solution for C 9i x,qU = 0 if P p (a + hr) = 0. Then, with the 
choice of ip and ip made above, this equation is equivalent to the following 

e v / h h 2m Cg >Xiqe -' t, /h(e-i'l , /h hr ) = - e -^ /h (h 2m /: g ^x,qa + h 2rn - 1 (Vp,X) g a). 

This will be solved by using Proposition 2.3. We find r £ H 1 (M) satisying 

IMIi^cdM) = 0 ( 1 ). 

The discussion of this section can be summarized in the following proposition. 

Proposition 3.1. Assume that ( M,g ) satisfies (9) and (10), and let m > 2 be an 
integer. Suppose that X is a W 1 ’ 00 vector field on M and q £ L°°(M). Let u> £ D 
such that (xi,u>) ^ M for all x\. If (r,6) are polar normal coordinates in (D,go) 
with center u>, then the equation 

Lg,X,qU = 0 in M 

has a solution of the form 

u = e~^ ( -‘ p+l ^\\g\~ 1 ^ 4 c 1 ^ 2 ao(xi : r, 6)b(6) + hr), 

where dao = ca\ for some a\ depending on (x\ , r) and satisfying da\ = 0 , b is 
smooth and the remainder term r £ such that ||r||jyi ^m) = 0(1). 

Remark 3.2. In fact, we need complex geometric optics solutions belonging to 
H 2m (M). Such solutions can be obtained in the following way. Extend X and q 
smoothly to K x Mq. By elliptic regularity, the complex geometric optics solutions 
constructed as above in M 0 will belong to H 2m (M). 

Remark 3.3. It is easy to check that if ao depends only on (x\ , r) and satisfies 
dao = 0, then the equation C g ,x,qU = 0 in M has a solution as in Proposition 3.1. 


4. Proof of Theorem 1.3 

Let (M, g) be an admissible manifold and let to > 2 be an integer. The first 
ingredient in the proof of Theorem 1.3 is a standard reduction to a larger compact 
manifold with boundary. 


Proposition 4.1. Let M, Mi be compact manifolds with boundary such that M CC 
Mi, and let m > 2 be an integer. Assume that Xi,X 2 are W 1,oc vector fields on 
M and qi,q 2 € L°°(M). Suppose that 


Xi=X 2 , qi = q 2 in Mi\M. 


If C™ x a = C^ x a , then Cfj- „ 


the Cauchy data for C g ^ Xj ,qj in Mi, j = 1,2. 


= r , , where Cfb 

(?,A2,<3'2 7 <?,A 


3 iQj 


denotes the set of 
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Proof. Let u £ H 2m (Mi) be a solution of C g> x 1 ,q 1 u = 0 in Mi. Since C^ Xl qi = 
C^ x 2 q 2 i there exists v £ H 2m [M ), solving C g X2 , q2 v = 0 in M, and satisfying 
jv = "fu in dM and jv = 7 u in dM. Setting 

f v in M, 

Vl = < 

( u in Mi \ M, 

we get Vi £ H 2m (Mi) and C g ,x 2 ,q 2 v 1 = 0 in Mi. Thus, C q % iqi C C™k 2 ,q 2 - Exactly 
the same way but in the other direction finishes the proof. □ 


The second ingredient is the derivation of the following integral identity based on 
the assumption that C™ Xiqi = 

Proposition 4.2. Let (M,g) be a compact Riemannian manifold with boundary, 
and let m > 2 be an integer. Assume that Xi,X 2 are W 1 ’ 00 vector fields on M and 
< 71,92 £ L°°(M). If C g , Xl , qi = C g ,x 2 ,q 2 , then 

/ [(Xi - X 2 , vVu) g + (91 - q 2 )uv} dVol g (x) = 0 , 

JM 

for any u,v £ H 2m (M ) satisfying C g ,-x 1 ,~ div g x^+ qi v = 0 and £ g ,x 2 ,q 2 n = 0 in 

M. 


Proof. We will use the following consequence of the Green’s formula, see [ 8 ], 

{Cg,Xi, qi u,v) L 2 ( M) = {u,C* Xi qi v) L 2 (M ) (11) 

for all u,v £ H 2m (M) such that ju = "/v = 0, where £* Xli?1 = £g,-x 1 ,-div a x 1 + qi - 


Now, let u, v £ H 2m (M) be such that C g -x x - div g A'l+qW = 0 and C g ^x 2 ,q 2 u = 0 in 
M. The hypothesis that M g ^Xi, qi = N g ,x 2 , q2 implies the existence of u £ H 2m (M) 
such that Cg,Xi , 9l u = 0 and 7 u = 7 u, f/u = 7 u. We have 

H-g,X x ,q x (u - u) = (Xi - X 2 )u + (91 - q 2 )u. 

Using (11), this implies the result. □ 


According to hypothesis, that X 1 = X 2 in (RxM 0 )\M lnt . We also extend 91 and q 2 
toRxM 0 by zero outside M lnt . Let, as in Section 3, ( D,g 0 ) CC (Mg nt , go) be simple 
such that ( M,g) CC (R x D lnt ,g) C (R x Mg 11 *, g). Let (Mi, < 7 ) be also admissible 
and simply connected such that ( M,g ) CC ( M[ nt ,g ) and (Mi,g) CC (M x D lnt ,g). 
According to Proposition 4.1, we know that C q I Xi = q 2 * s true. 

According to Proposition 4.2 the following integral identity holds for all u,v £ 
H 2m (Mi) satisfying C g ,x 2 ,q 2 u = 0 and C g _x u -dw g Xl +qiV = 0 in Mi, respectively: 



X 2 ,vVu) g + (91 


q 2 )uv] d Volg(x) = 0. 


( 12 ) 


The main idea of the proof of Theorem 1.3 is to use the integral identity (12) 
with u,v £ H 2m (M) being complex geometric optics solutions for the equations 
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C g ,x 2 ,q 2 u = 0 and Cg } -x 1 ,-di\„ Xi+qiV = 0 in Mi, respectively. We use Proposi¬ 
tion 3.1, Remark 3.2 and Remark 3.3 to choose solutions of the form 

u = e-^ xl+ir) (\g\- llA c 1/,2 e iX{xi+ir) b{6) + hr i), 

t, = es(* 1+ir )(| 5 |-i/4 c i/2 + /ir 2 ) 1 


where A G R and ||rj||#i ^ Ml ) = 0(1), j = 1, 2. Substituting these solutions in (12), 
multiplying the resulting equality by h and letting h —> 0 , we get 

lim / (Xi — X 2 , Xp) g uvd\o\ g (x) = 0, 

h ^° J Mi 

where p = x i + ir. Let us rewrite the integral in (x\,r,9) coordinates. Write 
X = X\ — X 2l and let be a 1-form dual to X. Let X xi and X\, denote the 
components of in the X\ and r coordinates. Then 


Mi. 




iX b r )e iX ( Xl+ir) b(6) dr dO dx\ = 0, 


(13) 


where M 1>X1 = {(r, 9) : {x\,r,9) G Mi}. Since X\ = X 2 in (R x M 0 )\M lnt , we may 
assume that the integral is over R xD. Taking aq-integral inside gives 


Define 



e lXxi (X»+iX°){ Xl ,r,6)dx 1 


dr d6 = 0 . 


f(x') = / e‘ All 4 1 (ii,a; , )(la:i, a(x') = ^ 

./R - 9 


e iAxi ^(s 1 ,iB , )da;i 


dad. 


Then / G W 1 , 00 (D) and a is a 1-form which is W 1,0 ° on D, and the integral identity 
above can be rewritten as 


[ [e Xr [f(7w,o(r))+ia('y w j(r))]drd9 = 0, 

Js n ~ 2 J 

where J w ,e is a geodesic in (D,go) issued from the point lo in the direction 9. For 
u! G dD , the integral above is related to the attenuated ray transform of function 
/ and 1-form ia in D with constant attenuation —A. Therefore, by varying the 
point ui in Proposition 3.1 on dD and using Proposition 5.1 in Section 5, for small 
enough A, we have / = —A p and a = —idp where p G W 1,0 °(D) and p\go = 0. The 
definition of a and analyticity of the Fourier transform imply that 

d k X\-djX ^ = 0 , j,k = 2,..., n. 

Also 

J e LXxi (djX\ — d\X'j){x\,x')dx\ =djf + i\atj =0, 

showing that dX b = 0 in Mi. Since Mi is simply connected, there is <j> G W 2, °°(Mi) 
such that 4>\dMi = 0 and X = V</>. 

Since X = X\ — X 2 in the neighborhood of the boundary dM\, we conclude that 
0 is a constant, say c G C, on dM\. Therefore, considering <f> — c, we may and will 
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assume that 0 = 0 on dM\. Since X\ — X 2 in (R x Mo) \ M lnt , we also may and 
shall assume that 4> is zero outside M\. In particular, </> is compactly supported. 
Next, we show that Xi = X 2 . For this, using Proposition 3.1 and Remark 3.2, 
consider 


u = e -i (x ' +xr \\g\- 1/4 c 1/2 e lX{xi+ir) b{d) + hn), 
v = e i(*i+ir)(| 9 |-V 4 c i/ 2 ao + hr2 ^ 


where ao satisfies dao = c. Such a 0 can be constructed using Cauchy’s integral 
formula in [ 6 ] as 

ao(xi,r,9) = ao(p,0) = — f C ^ dz A dz, for all 9 £ S n ~ 2 , 

2tt z- p 

where p = x\ + ir, B is a bounded domain in the upper half plane H C C such 
that the map B x S n ~ 2 -)Ix M 0 , (x\ ,r,9) i-a (xi, exp^(r0)) covers Mi and the 
boundary dB is piecewise smooth. Here and in what follows, u £ D such that 
w £ Mi in Proposition 3.1. 

Substituting these solutions and Xi — X 2 = V</> in (12), multiplying the resulting 
equality by h and letting h —> 0, we get 


lim / (V0, Vp) g uv d Vol s (x) = 0, 
h ~>° J Mi 

where p = x 1 + ir. Rewriting the integral in [x\,r,9) coordinates and taking x\- 
integral inside, we obtain 

2 / ([ ( dcj)a 0 e iX{xi+ir) b{e)dxidr ) d9 = 0 . 

Js n ~ 2 \J 0 J R / 

Since <f> is compactly supported, integrating by parts, in (xi,r), gives 


0 = - 


d<f> aoe iX ( Xl+ir ^b(6) dx\ dr] d9 


s n ~ 2 \J 0 J R 


</>da 0 e lHxi+ir) b(9)d Xl dr) d9 


<Pce iX{xi+ir) b(9)d Xl dr) d9. 


(14) 


Set 

$ A (r, 0 ) = [ (t>ce lXxi dx 1 , 

Jr 

i.e. is the Fourier transform of tpc in aq-variable. Then (14) can be written as 



e Xr ^x{ r yu,e{r))b(9)drd9 = 0. 


By varying the point u in Proposition 3.1 on dD and using [5, Lemma 5.1], for small 
enough A, we have = 0. Since 4>c is compactly supported, its Fourier transform 
is analytic. Therefore we obtain, (f> = 0 which shows that X\ = X 2 . 
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To show that qi = < 72 , consider (12) with X\ = X 2 which becomes 


/ (91 — q 2 )uv dVolg(x) = 0 

JMi 


(15) 


holds for all w,i> £ satisfying C g ,x 2 ,q 2 u = 0 and C g -x 1 -div g x 1 + qi v = 0 

in Mi, respectively. Use Proposition 3.1, Remark 3.2 and Remark 3.3 to choose 
solutions of the form 


u = e~i ( - Xl+ir \\g\~ 1 ^ 4 c 1 ^ 2 e iX( ' Xl+ir ' 1 b(6) + hn), 
v = e^ Xl+ir) {\g\~ 1/4 c 1/2 + hr 2 ), 


where AgK and ||rj||#i ( Ml ) = 0(1), j = 1,2. Substituting these solutions in (15) 
and letting h —>• 0, we get 


0 i\(x!+ir) 


(9i 


Mi, 


q 2 )c(x 1 ,r, 9)b(9) dr dd dx 1 = 0. 


Taking si-integral inside and varying b gives 



e lX ( Xl+lr ) (qi — q 2 )c{xi,r,9) dx\ dr d6 = 0. 


Set 


Qx(r,0) = [ (qi - q 2 )ce lXxi dx i, 

Jr 


i.e. Q\ is the Fourier transform of (q\ — q 2 )c in Xi-variable. Then, as in the case of 
d>A, one can show that Q\ = 0 for all A small enough. We have extended q\ and q 2 
to R x Mq by zero outside M, which implies that q± — q 2 is compactly supported. 
Hence, Q\ is analytic. This together with Q\ = 0 for all A small enough, allows us 
to conclude that q\ = q 2 - 


5. Attenuated ray transform 


The aim of this section is to prove the following proposition which was used in the 
proof of Theorem 1.3. We will closely follow the arguments in [5]. 


Proposition 5.1. Let ( D,go) be an (n — 1)-dimensional simple manifold, 
f £ L°°(D) and a be a 1-form which is L°° on D. Consider the integrals 




[/(7w,e(r)) + (Xk{luA r ))itfi{r)]e Xr b(9)drdd, 


Let 


where (r,0) are polar normal coordinates in (D,go) centered at some ui £ dD, and 
t(u>, 0) is the time when the geodesic r ha (r, 9) exits D. If |A| is sufficiently small, 
and if these integrals vanish for all cu £ dD and all b £ C°°(S n ~ 2 ), then there is 
p £ W 1,0 °(D) with p\qd = 0 such that f = —A p and a = dp. 


This is related to the injectivity of attenuated ray transform acting on function and 
1-form on D. Let us introduce some notions and facts; see [21] for more details. By 
SD we will denote its unit sphere bundle SD := {(&, v) £ TD : = 1}. On 
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the boundary of D, we consider the set of inward and outward unit vectors defined 
as 


d+SD = {(x, v) £ SD : x £ dD, (v, v(x)) go (x) < 0}; 
d-SD = {(x, v) £ SD : x £ dD, (v, v{x)) ga ( x) > 0}, 


where v is the unit outer normal to dD. The geodesics entering D can be param¬ 
eterized by d+SD. For any (x,v) £ SD the first non-negative exit time of the 
geodesic l j x ,v, with x = 7 Xi „(0), v = 7 Xi „(0), will be denoted as t(x,v). Simplic¬ 
ity assumption guarantees that t(x,v) is finite for all (x,v) £ SD. We also write 
4>t{x,v) = (nx,v(t),ix,v(t)) for the geodesic flow. 

We endow the unit sphere bundle SD with its usual Liouville (local product) mea¬ 
sure dS 2 " -3 , and endow the bundle d+SD with its standard measure d£ 2 " -4 . By 
dcr x we denote the measure on S X D. 

Let / be a function and a be a 1-form on D. The geodesic ray transform of / and 
a , with constant attenuation —A, is defined as 



In Propostion 5.1, if / and a were a continuous function and 1-form, respectively, 
one could choose b(9) to approximate a delta function at fixed angles 0 and obtain 
that 



for all to £ dD and all 9 £ S n ~ 2 . This would imply that 

T\[f,a](x,v) = 0 for all (x,v) £ d+SD. 

We will use the following result from [3, Theorem 7.1]. 

Proposition 5.2. Let ( D,go ) be a compact simple manifold with smooth boundary. 
There exists e > 0 such that the following assertion holds for a real number A with 
|A| < e: If f £ C°°(D) and a be a smooth 1-form on D, then T\ [/, a](x, v) = 0 
for all (x,v) £ d+SD implies the existence of p £ C°°(D) with p\qd = 0 such that 
f = —A p and a = dp. 

The previous argument together with the above theorem proves Proposition 5.1 for 
smooth / and a. However, this requires / and a to be C^-smooth in D and it is 
not obvious how to do this when / and a are L°° on D. We resolve this problem 
by using duality and the ellipticity of the normal operator TfT\. 

In the space of functions on d+SD define the inner product 



where dp,(x,v) = (v, v) go ^dY, 2n 4 . Denote the corresponding Hilbert space and 
the norm by L 2 ^(d+SD) and || • \\L^(d+SD ), respectively. We will also write 


h^{x,v) = h{(f_ T ( x _ v) {x,v)), (x,v) £ SD, 


for h £ C°°(d+SD). 
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If F is a notation for a function space (C k , L p , H k , etc.), then we will denote by 
T(D) the corresponding space of pairs [/, a] with / a function and a a 1-form on 
D. In particular, C 2 (D) is the space of square integrable pairs [/, a], and we endow 
this space with the inner product 


([/i«],[/'. a'])c*{D) = [ (ff'+(a,a') go )dVol go . 

Jd 

Lemma 5.3. If f £ C°°(D), a is a smooth 1-form on D and h £ C^°((d+SD) lnt ), 
then 

{T\[f,a],h) L 2( d+SD ' ) = ([/, a\,Tfh) C 2( D ), 
where Tfh is a pair defined as 

T\h(x) = f h 1 p(x,v)e~ XT< ' x ~ v ' l dcr x (v),f v k h^(x, v) e - XT ( x ~ v ') da x (v) . 

JS X D J S X D 

Proof. By Santalo formula (see [21] or [2]) 

ct],h) L 2( d+S D) 


L 


t(x,v) 


[filx,v{t)) + a k (j XtV (t))-y k }V (t)]e xt dt h{x, v) dp{x, v) 

d+SD J 0 

[/(*) + a k {x)v k ]h*(x, v)e~ Xr dH 2n ~ 3 (x, v) 


' SD 


1 [ f( x ) ( [ h^(x,v)e Xt( ' x ’ v) da x (v) \ dVo\ g (x) 

Jd \J s x d J 

+ f a k (x) ( f v k h^(x,v)e~ XT< ' x ~ v ' ) da x (v)\ dVol g (x). 
Jd \Js x d J 


This proves the statement. 


□ 


Proof of Proposition 5.1. First, we extend (D,go) to a slightly larger simple man¬ 
ifold and to extend both / and a by zero. Then / and a are still in L °°, and in 
particular in L p for all 1 < p < oo. In this way we can assume that both / and a 
are compactly supported in D lnt . 


We let b also depend on ui and change notations to write the assumption in the 
form 




e +a k ('Yx,v{t)) A / XjV {t)]b(x,v)dtd(T x (v) = 0 


for all x € dD and b £ Cpf((d+SD) lnt ). Let D be a compact submanifold of 
D with boundary such that ( D,go ) is also simple, D CC D lnt and supports of / 
and a are compact subsets of D lnt . Note that a is L°° on D (and in particular 
being in L 2 on D) implies that in particular 5a £ H~ 1 (D). Then we obtain the 
solenoidal decomposition a = a s + dp on D , where 5 go a s = 0 and p £ Hq(D) 
with —A go p = 5a. Here 5 go a = V^c^, where V 9o is the covariant derivative 
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corresponding to the metric go- Extend p to D by zero, so that p G Hq(D) with 
p = 0 in D \ D. An integration by parts shows that we have 

[ [ e~ Xt [f{lxAt)) + MlxA*)) + a k(lxA t ))'Yx, v ( t )}b(x,v) dtda x (v) = 0 

JS X D JO 

(16) 

for all x G dD and b € C^((d + SD) lnt ). Next, we make the choice b(x,v) = 
h(x,v)/j,(x,v) f° r h G Cg°((<9 + S'Z)) lrlt ) and integrate (16) over dD and get 

(T\[f + Xp,a s ],h) L 2^ d+SD ) = 0 . 

We are now in the same situation as in the proof of Lemma 5.3, and using the 
Santalo formula implies 

([/ + Ap, a s ],T x h)c 2 (D) — 0 

for all h G Cq°(( d + SD) lnt ). Note that the last integral is absolutely convergent 
because / G L°°(D) and a is 1-form which is L°° on D, and also the previous steps 
are justified by Fubini’s theorem. 

It remains to choose h = T x [tp, /?] for ip G Cfi°(D lnt ) and /3 being C^-smooth 1-form 
in D lnt , to obtain that 

([f + \p,a s ] 1 T*T x [<p,p])c 2 (D)=0. 

Since T X T X is self-adjoint, we have 

(T* x T x [f + Xp,a s },[p,/3])c 2 (D) =0 

for all ip G Cfi°{D int ) and for all C^-smooth 1-form /3 in D int . Therefore, T*T x [f + 
Ap,a s ] = 0. By [10, Proposition 1], T X T X is an elliptic pseudodifferential operator 
of order —1 in D mt . Here, ellipticity of T X T X is in the sense that whenever f,a' 
are in L 2 (D lnt ) and T x T x [f',a'] = 0 and S go a' = 0, then f',a' are smooth. Since 
f + Xp and a s were compactly supported in D lnt , this implies that / + Ap and a s 
are smooth and compactly supported in D lnt . Hence / + Ap and a s are smooth in 
D and compactly supported in D lnt . Now we can use the argument for smooth / 
and a given above, together with Proposition 5.2 to conclude that / = — Ap — Xip 
and a = a s + dp = dip + dp for some ip G C°°(D) with ip\dD = 0. To finish the 
proof, it remains to show that p G W 1 ’ 00 )!?). But this is clear from dp = a — a s 
and from a is L°° on D and a s is C°° on D. □ 

6. Boundary Determination and proof of Theorem 1.4 

In this part we show boundary determination of the vector field X. For the gener¬ 
ality of the statement we will assume the knowledge of the Cauchy data set C 9t x,q ■ 
It is easy to see that when 0 is not a Dirichlet eigenvalue of C g ,x,q, knowledge of 
the Cauchy data set is equivalent to knowledge of the Dirichlet-to-Neumann map 
Ng t x,q■ Moreover, we can determine not only the boundary values of X , but also 
the boundary values of q. This is the following proposition. 

Proposition 6.1. Let (M,g) be admissible, and let m > 2 be an integer. Suppose 
that X is a C°° vector field on M and q G C°°(M). Then the knowledge of the 
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Cauchy data set C 9j x,q determines the boundary values of X and the boundary 
values of q. 

To prove this proposition, it suffices to show that for any p £ dM, C^ x determines 
X{p) and q(p). In the following we will consider this local problem. 


Fix a point p £ dM and let (x',x n ) be the boundary normal coordinates near p, 
where x' = (aq,. .. , x„_i). In these coordinates dM corresponds to { x n = 0}, the 
vector field X becomes the differential operator X = X J , and the metric tensor 
can be written as 

g = g a pdx a (g> dx^ + dx n ® dx n . 

Here the in the following we use the convention that Greek indices run from 1 to 
n—1 and Roman indices from 1 to n. Denote Dj = then the Laplace-Beltrami 

operator in the boundary normal coordinates takes the form 

-A g = D 2 n + iE{x)D n + Q 2 (x, D x >) + Qi(x, D x <) (17) 


with E, Qi,Q 2 given by 

E(x) =\g aft d n g ap , (18) 

Q 2 {x,D x ,)=g af} D a Dp, (19) 

Qi{x,D x >) = - i(^g aP d a (\og\g\) + d a g a P)Dp. (20) 


Next we would like to write the 2m order equation 

£g,x, q u = (— A g ) m u + Xu + qu = 0 in M, m> 2 (21) 

as a second order system. To this end, introduce 

ui=u, U2 = (-A g )u, ..., = (~A g ) m_1 n 

and let U = (u ±,..., u rn ) T . By a standard reduction, (21) can be written as a 
system of equations in JJ: 

£ Aii ,a 12 ,a 0 U := (—A g <S)I+iAu(x, D x ')+iA 12 (x)D n +A 0 (x))U = 0 in M, ( 22 ) 
where / is the m x m identity matrix, An(x,D x i), Ai 2 (x) and Aq(x) are defined 

by 


{ 0 

0 . 

. 0\ 

, A 12 (x) := 

( 0 

0 . 

. 0\ 

0 

0 . 

. 0 

0 

0 . 

. 0 

V X a (x)D a 

0 . 

. o) 


l X n {x) 

0 . 

. oj 


Au{x, D x >) : 
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/ 0 -1 0 
0 0-1 


Ao (x) ■■= 


0 

V q(x) 


o 

o 


o 

o 


0 \ 
0 


-1 

0 


The associated Cauchy data set to the system (22) is 


Cau,a i2 ,a 0 := {(U\dM, 9 v U\om) ■ C-a 11 ,a 12 ,a 0 U ~ 0 in M, 


U G ( H 2 (M)) m }. 


It is easy to see that Ca ±1 ,a 12 ,a 0 and C 9t x,q are mutually determined, hence it suf¬ 
fices to show Ca 11} a 12 ,a 0 determines An, Ai 2 and Aq at p G dM. 


The following result gives a factorization of the operator Ca x1 ,a 12 ,a 0 - Similar tech¬ 
niques are employed in [3, 12, 16, 18]. 

Proposition 6.2. There is a matrix-valued, pseudodifferential operator B(x,D x >) 
of order 1 in x', depending smoothly on x n , such that 

£aii,a 12 ,a 0 = (D n <g>I+iE(x)®I+iA 12 (x)-iB(x, D x '))(D n ®I+iB(x, D x >)) (23) 

modulo a smoothing operator. Moreover, the principle symbol of the operator B(x, D x i) 
is (x,f')I. Here E{x) and Q 2 (x, D x >) are given by (18) and (19) respectively. 

Proof. Plug (17) into (22) we have 

£ah,a 12 ,a 0 ={D 2 n + iE{x)D n + Q 2 (x, D x >) + Qi(x, D x >)) ® I 
+ iAn(x,D X ') + iA 12 (x)D n + A 0 (x). 

Comparing this expression with (23) gives the following constrains on B(x,D x >) 
modulo a smoothing operator: 

B 2 (x,D x >) +i[D n ® I,B(x,D x ')} - E{x)B{x,D x ,) - A 12 (x)B(x,D x >) 
=Qi{x,D x /) ®I + Qi(x,D x ,) ®I + iA n (x,D x ,) + A 0 (x). (24) 

Let b(x, f) be the full symbol of B(x, D X '), then (24) implies on the level of symbols 
that 

’Y' —d x ,bD‘?,b+d n b-E{x)b-A 12 (x)b = Q 2 (x,£')I+Qi(x,£')I+iAn(x,£')+Ao(x). 

Z — J a\ s 
| a |>0 

(25) 

Let b ~ Sj<i bj where bj(x,(f ) is an m x m matrix with entries homogeneous of 
degree j in £' = (£ 1 ,... ,£ n _i). Collecting the terms homogeneous of degree 2 in 
(25) yields 

b 2 (x,a = Q2(x,ai, 

from which we can choose 


bi(x,?) = -y/Q 2 (x,?)I- 


(26) 
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Collecting the terms homogeneous of degree 1 in (25) yields 

Mi + Mo + ^2 + d n bi - E(x)bi - A 12 (x)bi = Qi(x,£')I + iAn(x,f'). 

|a|=l 

(27) 

Since &i (#,£') has been determined above, E(x) and Qi(x,£') are known from (18) 
(20), by some elementary linear algebra there exists a unique bo(x, £') satisfying this 
identity. Next collecting the terms of homogeneous of degree 0 in (25) implies 

b 2 0 + b 1 b. 1 + b. 1 b 1 + dZhd$bo+ J2 d Z h * d t ' h i + J2 l d >^ d P b i 

|o:| = l |a| —1 | Q£ [ ==2 

+d n b 0 - E(x)b 0 - A 12 (x)b 0 = A 0 {x ), (28) 

From which we can solve for b—i(x,£'). In general, the term bj(x , £') can be deter¬ 
mined by considering the terms homogeneous of degree j +1 in (25) . This completes 
the proof. □ 

Proof of Proposition 6.1. Using a similar argument as in [16, Proposition 1.2], we 
conclude that the Cauchy data set Ca 11 ,a 12 ,a 0 determines the operator B(x\ 0, D x >) 
modulo a smoothing operator. Consequently each bj\ Xn =o is determined, j < 1. It 
follows from (27) that the following expression is determined by the Cauchy data 
set Ca 1± ,Ai 2 ,Ao : 

-x n k=o VQ 2 (x', 0, e) + ix a \ Xn=0 £ a , ('er- 1 . 

Varying £' determines V"| Xn= o and X a \ Xn= o, a = 1,..., n — 1. Evaluating (28) on 
{x n = 0} shows that the Cauchy data set C i a 11 ,a 12 ,a 0 determines Ao\ Xrl -o , hence 
Q\x n =0‘ □ 

Proof of Theorem l.f. If 0 is not a Dirichlet eigenvalue of £ 9l Xi,gi and C g .x 2 ,q 2 i 
then N gt x liqi = X g ,x 2 ,q 2 implies C fh Xi, qi = C g ,x 2 ,q 2 - By Proposition 6.1 we con¬ 
clude that X\ = X 2 on dM. The result then follows from Theorem 1.3. □ 


7. Proof of Theorem 1.5 


We will follow the argument of [12, Theorem 1.3] and [3, Theorem 4]. Proceeding 
as in the proof of Theorem 1.3, 

We can derive the following integral identity (see (13)) 



(X h Xi +iX b r )e i ^ Xl+ir) b(8)drdddx 1 = 0. 


(29) 


This is similar to (13) but this time the integral is over M instead of 1 x M- ixi 
since we cannot extend the vector fields X\ and X 2 any more. Varying the smooth 
function b{9 ) leads to 



for all 6 G S n ~ 2 
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where Mg := {{x\ ,r) £ R 2 : (xi ,r,6) £ M} and p = x\ +ir. Define 

f{x') = f e lXxi Xx 1 (xi,x') dxi, a(x') = Y'' f f e lXxi X^(xi,x') dx\ \ dxT 

M j _2'^'® / 

Here we extend X as zero outside of M so that the integral in x\ can be over R. 
The above argument shows that 

j + ia{ A i u]fi {r))}dr = 0 for all 9 £ S n ~ 2 , 


where the r -integrals are integrals over geodesics 7^0 in n(M) C M 0 . Observe 
that a(x') is W 1,0 ° on 7 r(M) and f(x') is L°° on Under the assumption that 

go) is a simple (n — l)-dimensional manifold, we can apply Proposition 5.1 to 
D := 7 r(M) and conclude that for small enough A, we have / = —Ap and a = — idp 
where p £ W 1 , 00 (n(M )) and p\d-K(M) = 0. The definition of a and analyticity of the 
Fourier transform imply that 

d k X)-d 0 Xl= 0, j,k = 2,... ,n in M int . 


Also 

J e iXxi (.djX \ - 5iX})(a'i, x') dx 1 = d 3 f + iAaj = 0 in M int , 

showing that dX b = 0 in M. Since M is simply connected, there exists a function 
(j) such that V</> = X £ W 1 ’ oc (M). By [9, Theorem 4.5.12 and Theorem 3.1.7] we 
have (f) £ C ll, 1 (M). 


Next we need to show that (f> is constant on dM. In the case where we can extend 
A to be a compactly supported IU 1 ’ 00 vector field on a larger manifold, that (j) is 
constant on dM simply follows from the construction, but here we have to prove it. 
This is the content of the next proposition. 

Proposition 7.1. The function (j) is constant on the connnected boundary dM. 

Proof. Let us start by constructing more complex geometric optics solutions. From 
Proposition 3.1, Remark 3.2 and Remark 3.3 we can choose complex geometric 
optics solutions of the form 

u = e-i^-^XM-V^aoixur, 9)b{9) + h n ), 

v = e ^+ lr X\g\- 1 / A c 1 / 2 + hr 2 ), 

where dao = 0, A £ R and ||rj|| ff i ( (Mi) = 0(1), j = 1, 2. Note that in the previous 
construction we choose ao = e lX ( Xl+lr ^ but this time we need more ao’s. Substituting 
these solutions in ( 12 ), multiplying the resulting equality by h and letting h —> 0 , 
we get 
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where p = x\ + ir and X = X\ — X 2 . Recall that X := S7(f>. Insert the above 
complex geometric optics solutions yields 

I 5(j)ao(xi,r,0)b(6) dr dO dx\ = 0. 

JM 

Varying the smooth function b{6) leads to 

f dcfaodp A dp = 0 for all 0 £ S n ~ 2 . 

J M e 

Integrating by parts and using that dao = 0 gives 

[ ifiaodp = 0 for all 0 £ S n ~ 2 

JdMg 


(30) 


and for every ao with dao = 0 . 

On the other hand, noticing that in solving the eikonal equation (7), we may choose 
tp = x 1 but if) = —r. Then we can construct complex geometric optics solutions of 
the form 


u = e h 


(\g\ 1/4 c 1/2 ao(xi,r, 0)b{9) + hri), 


v = e ^- ir \\g\- l / A c 1/2 + hr 2 ), 
where dao =0, A Gl and ||ry ||#i ( Mj ) = 0(1), j = 1, 2. Here 


d = 


1 


d 


d 


dx 1 dr 


Using a similar argument as in the preceding paragraph we can derive 

<f>dodp = 0 for all 9 £ S n ~ 2 


J 

JdM e 


and for every do with dao = 0. In particular we can choose do = ao where ao solves 
dao = 0. Then taking complex conjugate gives 


(j)aodp = 0 for all 9 £ S' 


n —2 


(31) 


I dM e 


Im (p aodp = 0 


JdMg 

Combining (30) and (31) we see 

/ Re (j> aodp = 0, 

JdMg 

for all ao with dao = 0. Using the argument in [4, Section 5] implies that Re 4>\dM g = 
F\dM g for some non-vanishing holomorphic function F on Mg. Observing that Im F 
is a harmonic function in Mg and Ira F\gM e = 0, we conclude that F is real-valued 
and hence is constant on each connected component of dMg. Varying 9 shows 
that Reef) is constant along dM. Likewise we can show Iitk j> is also constant along 
dM. □ 

Proof of Theorem 1.5. Since <f> = c for some constant c along dM , replacing <f> by 
0 — c if necessary, we may assume <f> = 0 on dM. The rest part of the proof is the 
same as that of Theorem 1.3. □ 


22 


YERNAT M. ASSYLBEKOV AND YANG YANG 


Acknowledgement: The authors would like to thank Prof. Katya Krupchyk 
for her suggestions on an earlier version of this paper. The authors are also deeply 
grateful to Prof. Gunther Uhlmann for his generous support relating to this project. 

References 

[1] A. Calderon, On an inverse boundary value problem, Comput. Appl. Math. 25 (2006), 
133-138. 

[2] N. S. Dairbekov, G. P. Paternain, P. Stefanov, G. Uhlmann, The boundary rigidity problem 
in the presence of a magnetic field, Adv. Math. 216 (2007), 535-609. 

[3] D. Dos Santos Ferreira, C. Kenig, M. Salo, G. Uhlmann, Limiting Carleman weights and 
anisotropic inverse problems, Invent. Math. 178 (2009), no. 1, 119-171. 

[4] D. Dos Santos Ferreira, C. Kenig, J. Sjostrand, G. Uhlmann, Determining a magnetic 
Schrodinger operator from partial Cauchy data, Comm. Math. Phys., 271 (2007), 467- 
488. 

[5] D. Dos Santos Ferreira, C. Kenig, M. Salo, Determining an unbounded potential from 
Cauchy data in admissible geometries, Comm. PDE 38 (2013), no. 1, 50-68. 

[6] H. M. Farkas, I. Kra, Riemann Surfaces, Second Edition, Graduate Texts in Mathematics, 
71, Springer-Verlag, New-York, 1992. 

[7] F. Gazzola, H.-C. Grunau, G. Sweers, Polyharmonic boundary value problems, Springer- 
Verlag, Berlin, 2010. 

[8] G. Grubb, Distributions and operators, volume 252 of Graduate Texts in Mathematics. 
Springer, New York, 2009. 

[9] L. Hormander, The analysis of linear partial differential operators. I. Distribution theory 
and Fourier analysis , Classics in Mathematics. Springer-Verlag, Berlin, 2003. 

[10] S. Holman, P. Stefanov, The weighted Doppler transform, Inverse Probl. Imaging 4 (2010), 
111-130. 

[11] V. Isakov, Completeness of products of solutions and some inverse problems for PDE, J. 
Differential Equations 92 (1991), no. 2, 305-316. 

[12] K. Krupchyk, M. Lassas, G. Uhlmann, Inverse Boundary value Problems for the Perturbed 
Polyharmonic Operator, Transactions AMS, 366 (2014), 95-112. 

[13] K. Krupchyk, M. Lassas, G. Uhlmann, Determining a first order perturbation of the bihar¬ 
monic operator by partial boundary measurements, J. Funct. Anal., 262 (2012), 1781-1801. 

[14] C. Kenig, J. Sjostrand, G. Uhlmann, The Calderon problem with partial data, Ann. of 
Math. (2) 165 (2007), no. 2, 567-591. 

[15] K. Krupchyk, G. Uhlmann, Uniqueness in an inverse boundary problem for a magnetic 
Schrodinger operator with a bounded magnetic potential, Comm. Math. Phys., 327 (2014), 
993-1009. 

[16] J. Lee, G. Uhlmann, Determining anisotropic real-analtic conductivities by boundary mea¬ 
surements, Comm. Pure Appl. Math., 42 (1989), no. 8, 1097-1112. 

[17] R. Michel, Sur la rigidite imposee par la longueur des geodesiques, Invent. Math., 65 
(1981), 71-83. 

[18] G. Nakamura, Z. Sun, G. Uhlmann, Global identifiability for an inverse problem for the 
Schrodinger equation in a magnetic field, Math. Ann. 303 (1995), no. 3, 377-388. 

[19] M. Salo, Inverse problems for nonsmooth first order perturbations of the Laplacian, Ann. 
Acad. Sci. Fenn. Math. Diss. 139 (2004). 

[20] M. Salo, G. Uhlmann, The attenuated ray transform on simple surfaces, J. Diff. Geom. 88 
(2011), no. 1, 161-187. 

[21] V. A. Sharafutdinov, Integral geometry of tensor fields, Inverse and Ill-Posed Problems 
Series, VSP, Utrecht, 1994. 

[22] Z. Sun, An inverse boundary value problem for Schrodinger operators with vector poten¬ 
tials, Trans. Amer. Math. Soc. 338 (1993), no. 2, 953-969. 

[23] J. Sylvester, G. Uhlmann, A global uniqueness theorem for an inverse boundary value 
problem, Ann. of Math. (2) 125 (1987), no. 1, 153—169. 


POLYHARMONIC OPERATOR ON ADMISSIBLE MANIFOLDS 


23 


[24] C. Tolmasky, Exponentially growing solutions for nonsmooth first-order perturbations of 
the Laplacian , SIAM J. Math. Anal. 29 (1998), no. 1, 116-133. 

[25] Y. Yang, Determining the first order perturbation of a bi-harmonic operator on bounded 
and unbounded domains from partial data , J. Differ. Equations 257 (2014), 3607-3639. 

Department of Mathematics, University of Washington, Seattle, WA 98195-4350, USA 
E-mail address : y_assylbekov@yahoo.com 

Department of Mathematics, Purdue University, West Lafayette, IN 47907, USA 
E-mail address : yang926@purdue.edu 


